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ABSTRACT 

■ We investigate properties of r-modes characterized by regular eigenvalue problem in slowly 

I rotating relativistic polytropes. Our numerical results suggest that discrete r-mode solutions 

for the regular eigenvalue problem exist only for restricted polytropic models. In particular the 
. r-mode associated with I — m — 2, which is considered to be the most important for gravitational 

^! I radiation driven instability, do not have a discrete mode as solutions of the regular eigenvalue 

problem for polytropes having the polytropic index N > 1.18 even in the post-Newtonian order. 

, Furthermore for a = 1 polytrope, which is employed as a typical neutron star model, discrete 

I r-mode solutions for regular eigenvalue problem do not exist for stars whose relativistic factor 

' M/R is larger than about 0.1. Here M and R are stellar mass and stellar radius, respectively. 

> ■ 

^ ; 

, Subject headings: relativity — stars: neutron — stars: oscillations — stars: rotation 

o : 
^ , 

. 1. Introduction 

Since the discovery of the gravitational radiation driven instability of the r-modes by Andersson (1998) 
^ I and Friedman & Morsink (1998), a large number of studies on the properties of r-modes and inertial modes 

of rotating stars have been done to prove their possible importance in astrophysics (for a review see, e.g., 
^ , Friedman & Lockitch 1999; Andersson & Kokkotas 2000). Although early investigations were mainly 

' applied to young hot neutron stars (e.g., Lindblom, Owen, & Morsink 1998; Owen et al. 1998), recently 

j_j ' some authors have also discussed possible roles of the r-mode instability in old and cool neutron stars 

I with a solid crust and a magnetic field, such as those found in low mass X-ray binaries (LMXBs) (e.g., 

Andersson, Kokkotas, & Stergioulas 1999; Bildsten & Ushomirsky 2000; Rezzolla, Lamb, & Shapiro 2000; 
Yoshida & Lee 2001). 

Although our understandings of r-modes have been largely improved by recent energetic investigations 
as mentioned above, most studies concerning r-modes have been done in the framework of Newtonian 
gravity. As for r-modes in the framework of general relativity, Kojima (1998) found possible existence of 
continuous spectrum for pure r-mode in relativistic stars. Here throughout this paper we call the rotational 
mode whose eigenfunctions are composed of only one axial parity component in the non-rotating star 
limit "pure r-mode". Beyer & Kokkotas (1999) generally verified the existence of a continuous spectrum. 
Kojima's formalism was developed to that with high order rotational effects by Kojima & Hosonuma 
(1999; 2000). Very recently, Lockitch, Andersson, & Friedman (2001) have shown that pure r-modes in a 
relativistic star can exist only in non-isentropic stars, whose specific entropy distribution of a star is not 
constant. As for rotation induced modes in isentropic stars, they have shown that modes similar to inertial 
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modes, namely hybrid modes, in a Newtonian rotating star (see, e.g., Lockitch & Friedman 1999; Yoshida & 
Lee 2000), whose eigenfunctions are composed of both polar and axial parity components, exist in uniform 

density stars (sec also Lockitch 1999). 

In the study by Lockitch et al. (2001), they derived the basic equation for the pure r-mode in relativistic 
non-isentropic stars. Their equation is the same as that obtained by Kojima (1998). As mentioned by 
Kojima (1998), this equation has a continuous part of frequency as its solution when the system of equations 
is that of the singular eigenvalue problem. Kojima (1998) also pointed out that the system of the equation 
becomes that of the regular eigenvalue problem when the frequency domain is in some range. For such a 
regular eigenvalue problem, discrete eigenfrequencies may exist for the same equation as well as continuous 
ones. Lockitch et al. (2001) have studied the relativistic r-modc charactcirized by the regular eigenvalue 
problem. In their study discrete r-mode solutions for the regular eigenvalue problem in uniform density 
stars were obtained. It is reasonable to consider that such discrete r-mode solutions correspond to those of 
Newtonian r-mode. However their study was restricted to the case of uniform density stars. As mentioned 
in the next section, the frequency range for the regular eigenvalue problem is confined to the narrow region 
and depends on the structure of equilibrium stars. Furthermore, needless to say, relativistic effects on 
stellar oscillations in a neutron star are not negligible. Thus, it is interesting and important to answer 
to the question of whether r-mode solutions for the regular eigenvalue problem exist in all non-isentropic 
relativistic stars. In this paper, therefore, we study the r-mode characterized by regular eigenvalue problem 
for large class of polytropic and non-isentropic relativistic stellar models. The plan of this paper is as 
follows. We start, in §2, with brief description of our method of solution. In §3, we describe the properties of 
the r-mode characterized by the regular eigenvalue problem in both the relativistic and the post-Newtonian 
stellar models. §4 is devoted for discussions and conclusions. Throughout this paper we will use units in 
which c = G = 1, where c and G denote the velocity of light and the gravitational constant, respectively. 



2. Method of solutions 

We consider slowly rotating relativistic stars with a uniform angular velocity il, where we take account 
of the first order rotational effect in fl. The geometry around the equilibrium stars can be described by the 
following line element (see, e.g. Thorne 1971): 

rfs^ = -e'^'^^^^dt'' + e^^^^Ur'' + r^d0^ + sin^ edip^ - 2w(r) r^ sin^ Odtdip . (1) 

Once an equation of state, p = p{e, s) for a perfect fluid is given, we can easily obtain the structure of 
the slowly rotating stars, where p, e, and s denote pressure, mass-energy density and specific entropy, 
respectively. Throughout this paper, the polytropic equation of state is assumed: 

p = Ke^+7r, (2) 

where N and K are constants. 

We consider oscillation modes in rotating relativistic stars such that the eigenfunctions are stationary 
and are composed of only one axial parity component in the limit 17 ^ 0. Subclass of those modes should be 
a relativistic counterpart of r-modes, which are able to oscillate in all slowly rotating Newtonian fluid stars. 
According to Lockitch et al. (2001), such modes are allowed to exist only if the star has non-isentropic 
structures. Therefore, we assume stars to be non-isentropic, although the effects due to deviation from 
isentropic structure on oscillation modes do not appear in the first order in O. According to the formalism 
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by Lockitch et al. (2001) (see, also, Kojima 1998), let us write down the pulsation equation for relativistic 
r- modes with accuracy up to the first order in fl. The metric perturbation 6gaf3 and the Eulerian changes 
of the fluid velocity Su" that do not vanish in the limit ^ are given as 

iSgte,Sgt^) = iho,i{r) f- ^^^'"^f ' , sin 9deYirn{0, ^)) e'^' , (3) 



where Yi„i{0,ip) are usual spherical harmonic functions, and a denotes oscillation frequency measured 
in the inertial frame. All other perturbed quantities become higher order in Q. Note that this form of 
eigenfunctions is the same as that for zero-frequency modes in a spherical non-isentropic star, because 
r-modes become zero-frequency ones in the limit fl — *■ 0.(Thorne & Campolattaro 1967) The metric 
perturbation ho^i obeys a second order ordinary differential equation, 
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Dim{r\a) 

+ 167r(p + e)/io,i = , (5) 



where 



DUr-a) = l- , 6 

L{L + l)cr 

Here, we have introduced the effective rotation angular velocity of fluid, 

u) — Q. — uj , (7) 

and the corotating oscillation frequency, 

a = cr + mn . (8) 

The velocity perturbation of fluids Ui is determined from the function h^j through the following algebraic 
relation, 

2muj 



/(t -I- l)a 

Equations (^ and (|^) are our basic equations, which were derived by Kojima (1998) in his early work. Note 
that equations (^ and (^) lose the meaning for the I = case, because there are no axial modes with I = 0. 

Because equations (^) are second order ordinary differential equations, two boundary conditions are 
required to determine solutions uniquely. For physically acceptable solutions, the function /iq,; must vanish 
both at the center and at the spatial infinity. Thus, boundary conditions are given as 



dhp,, 

dr 
dhp^i 

dr 



- {l + l)hoi^O asr^O, (10) 

dr 

I ho^i = as r — > 00. (11) 



When those boundary conditions are imposed, our basic equations are solved as an eigenvalue problem with 
an eigenvalue a. As shown by Kojima (1998), when Di,n{r; ct) = is satisfied inside the star the system of 
our basic equations becomes a singular eigenvalue problem, because the function Dim{r\a) is proportional 
to the coefficient of the second derivative of function /iq,;. On the other hand, we can treat the equation 
as a regular eigenvalue problem if Dim{r]a) = is not satisfied inside the star, because the last term of 



-4- 



equation (|^) vanishes outside the star as weU as the function Di,n(r; a). Due to this singular property of the 
basic equation, Kojima also concluded continuous spectrum may exist in relativistic rotating stars. Beyer 
& Kokkotas (1999) gave a mathematically rigorous proof for the existence of a continuous part within the 
spectrum. Recently Lockitch et al. (2001) have showed that the equality Dim{r;a-) — must be satisfied 
at some spatial point for the equations to have non-trivial solutions hoj. Because w is a monotonically 
increasing positive function of r, according to the considerations of Lockitch et al. (2001), in order for the 
equation (^ to be treated as a regular eigenvalue problem an frequency a must be in a region, 

2muj{R) _ 2mLd{oo) 2mil 

1(1 + 1) 1(1 + 1) "1(7+1)' ^ ' 

where R is stellar radius in the coordinate (|l|). Note that is well-known Newtonian r-mode frequency. 

On the other hand, when a satisfies an inequality, 

2m^(0) . 2mu^(R) 

7(7TT) -"-7(7TTy' ^^'^ 

we must treat equation (|^) as singular eigenvalue problem. Note that since u) becomes fl in Newtonian 
limit, we can easily see from inequalities (^2|) and ( p^ that frequency in a range for existing non-trivial 
solution of equation (^ is getting closer to Newtonian r-mode frequency, as M/R 0. Lockitch et 

al. (2001) also found numerically discrete r-mode solutions with frequencies in range ( [l^ ) in uniform density 
relativistic stars. In this paper we will extend the study on the discrete r-mode solutions by Lockitch et al. 
(2001) to those for polytropic equation of states. 

Next let us consider a post-Newtonian approximation of relativistic r-mode oscillation, by using a 
similar method to that of Lockitch et al. (2001). In this approximation, physical quantities of an equilibrium 
state can be written as 

'-'-^'-¥^00)')^ 

^'A^--'^^^'-¥)^o0)y 

where is a Newtonian Lane-Emden function (see, e.g., Kippenhahn & Weigert 1990) and M denotes total 
mass of stars. Here eg and Pq are constants determined from a solution of Lane-Emden equation. Function 
Uq is given as solution of a differential equation, 

^^ + ^^ + 16.|Le- = 0, (19) 
dr^ r ar /t^ 

under the following boundary conditions: 

^=0, u.o(R) + ^^^=0. (20) 
dr 3 dr 
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In the Newtonian limit, — > 0, eigensolutions of equation (||) should become those of r-mode 
oscillations in Newtonian stars. Therefore, eigenvalue a and eigenfunction can be expanded in power series 



oi^ as 



2mVL 
l{l + l) 



,2M , ^ f2M\' 



and 



Ui = ul°^ + o 



(2M\ 



\ R J 

Substituting equations and ( pT| ) into an equation (||), we obtain basic equations for r-mode 

oscillations in the post-Newtonian order as follows: 



(21) 
(22) 

(23) 



(ki +ujo) 



^2 "-0,1 V) 



16.|Le^4°)(r) = 0, 



(24) 



and 



(0) 



,(0) 
'■0,1 



Kl + UJQ 

These equations have the same properties as those of relativistic equations (^. Boundary conditions 
required for physically acceptable solutions are given as follows: 



(25) 



dh, 



(0) 
0,1 



{I + 1) h'^] =0 as r ^ , 



dr 



and 



dh, 



(0) 

0,; 



dr 



lh^°} =0 at r = i? . 



From inequalities ( p^ ) and (13), frequency range for regular and singular eigenvalue problem in the 
post-Newtonian order equation are given as 



and 



respectively. 



-UJq{R) < Kl < 0, 
-WO(0) < All < -LUo{R), 



(26) 
(27) 

(28) 
(29) 



Because our basic equation for both relativistic and post-Newtonian r-modes become second order 
linear ordinary differential equations with two boundary conditions, those are numerically solved as an 
eigenvalue problem, by using a Henyey type relaxation method (e.g., Unno et al. 1989). 



3. Numerical results 

3.1. Relativistic study 

As mentioned in the last section we should distinguish two cases in treating equation (^). One is 
a regular eigenvalue problem. The other is a singular eigenvalue problem. In this study we consider 
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only r-mode solutions of a regular eigenvalue problem, because discrete eigenfrequency may be obtained. 
Furthermore such discrete modes are considered to correspond to Newtonian r-modes, because such modes 
have the well-defined regular eigenfunctions. Therefore we restrict our attention to eigenfrequency in 
the range ([l^). For simplicity, henceforth, we will naively call a r-mode solution of a regular eigenvalue 
problem a discrete r-mode unless we give a notice. As mentioned in Beyer & Kokkotas (1999) we cannot 
exclude the possibility that isolated eigenvalues may exist within a frequency range for a singular eigenvalue 
problem. Thus note that such an isolated eigenvalue does not belong to "discrete r-mode" according to the 
terminology of this paper. 

We compute frequencies of a discrete r-mode for several polytropic stellar models. In present study, 
only the fundamental r-modes, whose eigenfunction Ui has no node in radial direction except at stellar 
center, are obtained. This is consistent with that in Newtonian case, where overtone r-modes appear as 
high order effects in angular velocity (see, e.g. Provost, Berthomieu & Rocca 1981, Saio 1982). In Figure 
1, scaled eigenfrequencies n = a/fl of discrete r-modes are given as functions of M / R for modes with ni = I. 
Note that for modes associated with /, we can directly obtain an eigenfrequency of a mode with any value of 
m from that with m = I. Therefore only results for modes associated with m — I will be shown throughout 
this paper. Eigenfrequencies for stars with four different polytropic indices, A'^ 0, 0.5, 0.75, and 1, are 
shown in panels in Figure 1, respectively. In each panel, frequency curves for the modes with ^ = 2, 3, and 
4 are depicted along a relativistic factor M / R. From the figure it is found that frequencies n monotonically 
decrease as a relativistic factor MjR increases. This behavior is caused by dragging of inertial frame due 
to the stellar rotation. We find that some frequency curves in Figure 1 are terminated at some value of 
MjR beyond which equilibrium states can still exist. Here, the maximum values of MjR for polytropic 
equilibrium stars having iV = 0, 0.5, 0.75, and 1.0 are given by 4/9, 0.385, 0.349, and 0.312, respectively. 
It is also found that length of those curves tend to be short as either a polytrope index N or an angular 
quantum number / increases. As a matter of fact, we are not able to obtain eigenfrequencies beyond them 
by using our numerical method. 

In order to trace the cause that missing frequencies in Figure 1 cannot be obtained, let us see behavior 
of values of ;™+i)o ^ ^ function of MjR. The value of gives a lower bound of the eigenfrequency k 

allowed for the regular eigenvalue problem of a r-mode in the considered star. Therefore an eigenfrequency 
K for the regular eigenvalue problem must be larger than . In Figure 2, values of and k for 

I = m = 2 r-mode of A^ = 0.75 polytropic stars are given as function of M/R. In this figure it is found that 
the curve of k approaches that of twice, that is, near the Newtonian limit, and near the terminal 

point of the frequency curve. The former is clear because the frequency range (|l2|) shrinks to the point of 
the Newtonian r-mode frequency. The latter is considered to be the cause of our problem. From the figure 
it is found that frequency n of r-mode solutions for the regular eigenvalue problem becomes the same value 
as that of ^T+i)^^ some critical value of M/R, say {M/ R)crit, within the numerical accuracy. In other 
words, K nearly becomes the frequency for the singular eigenvalue problem as the value of M/R becomes 
(M/ R)crit- Thus we can not obtain those missing modes. Behaviors of 'ff^^^^ and n as functions of M/R 
for other cases except for the N = case are similar to those in Figure 2. As for the A^ = case, the curve 
of K does not closely approach that of even when the value of M/R is nearly 4/9. From those facts 

and continuity of physical quantities we can expect that disappearance of r-mode solutions with frequency 
in a range ( p^ ) is not due to numerical artifact. Thus it is reasonable to claim that r-mode solutions for 
regular eigenvalue problem disappear in stars having a values of M/R larger than {M / R)crit- 

Critical values of M/R are shown as functions of polytropic index A^ in Figure 3. Curves of critical 
values {M/R)crit for modes with / = 2, 3, and 4 are given in the figure. Note that values of {M/R)crit shown 
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here are numerically approximated ones. These only give lower limit of accurate values of {M / R)crit- In 
Figure 3, stars whose parameters are in regions under curves of (M/ R)crit have discrete r-mode spectrum. 
Results of our numerical analysis also suggest that for the N = case the critical value oi M/R is nearly 
equal to 4/9. Indeed, discrete r-mode solutions with frequency in a range ( p^ may exist in all uniform 
density stars independently of the value of M/R. This is consistent with the result by Lockitch et al. 
(2001). The figure also suggests that even in nearly Newtonian stars, whose value of relativistic factor M/R 
is sufficiently small, r-mode solutions of a regular eigenvalue problem disappear if a value of N is over some 
critical value, say Ncrit- For such nearly Newtonian stars, however, our numerical method tends to lose 
an accuracy. In the next subsection, to avoid this numerical difficulty, we will study r-mode oscillations of 
nearly Newtonian stars by using post-Newtonian formalism. 



3.2. Post-Newtonian study 

In post-Newtonian analysis of r-modes, the post-Newtonian correction ki for the eigenfrequency is 
employed as the eigenvalue instead of a scaled eigenfrequency k. The value of ki does not depend on a 
relativistic factor M/R by definition but both on polytropic index TV and on angular eigenvalue /. In 
similar way to the case of the relativistic study we consider only the case that our basic equations become 
the regular eigenvalue problem. Thus we restrict eigensolutions to those with frequency in a range (p8|). 
In Figure 4 the post-Newtonian correction ki are plotted against polytropic index TV for three different 
eigenmodes associated with Z = 2, 3 and 4. Our numerical results of ki for a = polytropic model are 
in good agreement with those by Lockitch et al. (2001). We also cannot obtain r-mode solutions for the 
regular eigenvalue problem in stars beyond some limits of polytropic index N , as shown in the figure. Thus 
as in the relativistic case, examining behaviors of values of —ujq(R) appearing in the frequency inequality 
( p8| ) as a function of N, we confirm that disappearance of r-mode solutions with frequency in a range ( |2^ ) 
is not due to numerical artifact. The values of —uJo{R) are shown as function of iV in Figure 4. It is found 
that all frequency curves for ki in Figure 4 nearly intersect with the curve for —u!o{R) at some values of iV, 
namely N = Ncrit- In the same way as the relativistic study, thus, we can conclude that the absence of the 
post-Newtonian r-mode solutions for the regular eigenvalue problem in stars having iV > Ncrit is not due to 
numerical artifact. From this we expect that r-mode solutions with a frequency in a range ( p8|) disappear 
if a star has polytrope index N larger than a critical value Ncrit ■ Values of Ncrit obtained from numerical 
analysis are shown in Figure 4 near the terminal points of ki curves. Our numerical results of Ncrit in 
the post-Newtonian order are consistent with those in relativistic analysis, (compare with Figure 3) From 
the figure, for instance, we can expect that the r-mode solution with I — 2 for regular eigenvalue problem 
disappears in polytope having iV > 1.18 even in the post-Newtonian order. 



4. Discussion and Conclusion 

In this paper, we have investigated the properties of r-modes characterized by regular eigenvalue 
problem in a slowly rotating relativistic polytrope. Our numerical results suggest that discrete r-mode 
solutions with frequency in a range (|l^) disappear for some polytropic models. Although we do not 
have mathematically rigorous proof for this suggestion, from numerical results and continuity of physical 
quantities we can make conjectures on the relativistic r-mode as follows: the r-mode associated with 
I ^ m = 2, which is considered to be the most important for gravitational radiation driven instability, 
do not exist as a solution for regular eigenvalue problem for polytropes having iV > 1.18 even in the 
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post-Newtonian order. As for a N — I polytropic model, which is employed as a typical neutron star model, 
discrete r-mode solutions with frequency in a range do not exist for stars whose relativistic factor 
M/R is larger than about 0.1. If those conjectures are true, most of typical neutron stars do not have 



discrete r-mode solution with frequency in a range (12), because a typical relativistic factor for a neutron 



star is considered to be about 0.1-0.2. This means that in actual neutron stars r-mode characterized by 
regular eigenvalue problem cannot play any important role, if relativistic r-mode oscillations are governed 
by equation (^. If this is the case, we must study the pure r-mode with frequency in a range ( [Tsl ) to have 
clear understanding of relativistic r-mode in neutron stars. According to previous studies, continuous part 
of spectrum and singular eigenfunctions appear for such a frequency range, although we cannot exclude the 
possibility that isolated eigenvalues may exist within the frequency range. 

In this study, we employ the formalism developed by Lockitch et al. (2001) to investigate r- modes in 
relativistic polytropes. Although their formalism is correct, the system of equation (^) as regular eigenvalue 
problem leads to some unnatural numerical results. For instance, I — 2 r-mode solutions of regular 
eigenvalue problem in iV = 1 polytropes having AI/R > 0.1 disappear due to general relativistic effect, 
although r-mode solutions for regular eigenvalue problem in Newtonian stars exist for all slowly rotating 
polytropic stars with constant angular velocity. Lockitch et al. (2001) claimed that r-mode solutions with 
frequency in a range (|l2|) are counterpart of Newtonian r-modes, because they considered that the true 
relativistic analogue to a Newtonian r-mode should be a distinct mode with a well-defined eigenfunction. 
In fact, such a r-mode solution may exist in all A'^ = polytrope, as shown numerically in Lockitch et al. 
(2001) and this study. However, for some N ^ polytropes we numerically found that r-mode solutions 
with frequency in a range (|l^ ) do not appear when AI/R > {AI/R)crit or A > Ncrit- We think that 
this is actual behavior of the mode due to the influence of the effective differential rotation caused by 
general relativity. The reason of this is considered as follows: In modal analysis in Newtonian rotating 
stars, sufficient differential rotation of an equilibrium star can give a corotation point, where corotation 
frequency a vanishes. Because of existence of this corotation point the system may become singular 
eigenvalue problem, (see, e.g. Karino, Yoshida & Eriguchi 2000) On the other hand, our results show that 
disappearance of discrete r-mode occurs in large values of the relativistic factor or the polytropic index of 
the equilibrium star. In such equilibrium stars, as a matter of fact, differential rotation effect in relativistic 
rotating stars becomes strong. 

However we cannot exclude the possibility that pure r-modes in a relativistic rotating star with the 
non-zero rotation frequency are not simply governed by equation (H) , because strictly speaking the equation 
(^ only determines zero frequency modes in a spherical non-rotating star. In other worlds, in equation (^) 
the gravitational radiation reaction is not taken into account at all, although the gravitational radiation 
must be emitted from r-mode oscillations in rotating stars with non-zero rotation velocity. As strongly 
suggested by Lockitch et al. (2001) (see, also Beyer & Kokkotas 1999), when the eigenfrequency becomes 
complex number, which expresses the damping of the oscillation due to the energy dissipation such as the 
gravitational radiation emission, the existence of the continuous part of spectrum in the eigenfrequency 
of equation (|^) might be avoided. If this is true, there is a possibility for a distinct mode solution with a 
well-defined eigenfunction from singular solutions with frequency in a range ( p^ ) to be picked up by the 
information due to the gravitational radiation reaction effect. Such distinct modes may make a continuous 
sequence that starts with r-mode solutions with frequency in a range (^2|). We think this one of the most 
plausible scenarios, although this is still speculation. In order to have clear understanding of the r-mode in 
relativistic stars, the further progress of investigations of the rotational mode is needed. 
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Fig. 1. — Scaled frequencies k = a/Q of the r-modes in polytropic models with N = 0, 0.5, 0.75, and 1 are 
plotted as functions of M/R in panels, respectively. In each panel, the frequencies for modes with I = 2, 3, 
and 4 are given. The labels indicating their polytropic index N are attached in corresponding panels. The 
eigenvalue of the r-modes that is described above but not given in the figure, for instance I = 3 mode in 
A'^ = 1 polytrope, are not obtained by using our numerical method. Vertical dotted lines show the maximum 
values of M/R for equilibrium states: M/R = 0.444 for N = 0, M/R = 0.385 for N = 0.5, M/R = 0.349 for 
N = 0.75, and M/R = 0.312 for N = 1.0. 
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Fig. 3. — The critical values of relativistic factor M/ R are plotted against the polytropic index N for 
discrete r-modes with / = 2, 3, and 4. The values of {M/R)crit obtained from actual numerical calculations 
are indicated by the filled circles. 
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Fig. 4. — The post-Newtonian correction ki for discrete r-mode frequencies with / = 2, 3, and 4 are shown 
as functions of the polytropic index N. The labels indicating the critical values of N are attached near the 
terminal points for each frequency curves. The value of —uio{R) is also plotted as a function of N. 



